arXiv:1505.05297v2 [gr-qc] 13 Oct 2015 


Gravitational field equations near an arbitrary null surface 
expressed as a thermodynamic identity 

Sumanta Chakraborty“,* Krishnamohan Parattu“,^ 

and T. Padmanabhan“^ 

“lUCAA, Post Bag 4, Ganeshkhind, 

Pune University Campus, Pune 411 007, India 

October 14, 2015 


Abstract 

Previous work has demonstrated that the gravitational field equations in all Lanczos- 
Lovelock models imply a thermodynamic identity TiIaS = 5\E-\-P5\V (where the variations 
are interpreted as changes due to virtual displacement along the affine parameter A) in the 
near-horizon limit in static spacetimes. Here we generalize this result to any arbitrary 
null surface in an arbitrary spacetime and show that certain components of the Einstein’s 
equations can be expressed in the form of the above thermodynamic identity. We also obtain 
an explicit expression for the thermodynamic energy associated with the null surface. Under 
appropriate limits, our expressions reduce to those previously derived in the literature. The 
components of the field equations used in obtaining the current result are orthogonal to the 
components used previously to obtain another related result, viz. that some components of 
the field equations reduce to a Navier-Stokes equation on any null surface, in any spacetime. 
We also describe the structure of Einstein’s equations near a null surface in terms of three 
well-defined projections and show how the different results complement each other. 


1 Introduction 

Horizons in general (and black holes in particular) possess thermodynamic attributes like entropy 
[1,2] and temperature [3-6]. These features, which are known to transcend Einstein’s gravity, 
are believed to stem from some deep connection between gravitational dynamics and horizon 
thermodynamics. In recent years, it has been shown that in general relativity as well as in a wider 
class of gravitational theories, the field equations near a horizon imply a thermodynamic identity 
T5\S = S\E + PS\V [7-30] where the symbols have the usual meanings and the variations are 
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interpreted as changes due to virtual displacement along the afhne parameter A. This result — 
originally obtained for general relativity [7] — has been generalized to all static spacetimes with 
horizon in the Lanczos-Lovelock theories of gravity [13,15,31]. That is, the result is known to 
hold for actual horizons — rather than to generic null surfaces — and requires the assumption 
of a static spacetime. 

It may therefore appear that this connection — between the field equations and a thermo¬ 
dynamic identity — is an exotic phenomenon that occurs only in specihc solutions containing 
horizons. But this illusion is broken when we realize that a generic null surface through any event 
in spacetime can act as a local Rindler horizon for some observer [9,32]. This fact allows one to 
introduce observer-dependent thermodynamic variables around any event in spacetime and rein¬ 
terpret the gravitational field equations near any null surface in a thermodynamic language. One 
can then ‘derive’ field equations in the case of Einstein’s theory [32] from the Clausius relation 
applied to a null surface, if one assumes further that (a) the entropy density is one quarter of the 
transverse area and, more importantly, (b) the quadratic terms in the Raychaudhuri equation — 
involving the squares of shear and expansion — can be set to zero. (This approach based on the 
Raychaudhuri equation, however, could not be generalized in a simple manner to more general 
class of theories.) A clearer connection between null surface thermodynamics and gravitational 
dynamics emerged from the fact that, gravitational held equations reduce to the Navier-Stokes 
equations of huid dynamics in any spacetime when projected on an arbitrary null surface [33,34], 
thereby generalizing previous results for black hole horizons [35]. 

In the light of these results, it is natural to ask whether gravitational dynamics is the long 
wavelength thermodynamic limit of the dynamics of some unknown microscopic degrees of free¬ 
dom. The conceptual framework that attempts to interpret gravitational dynamics as emergent 
from the dynamics of unknown microscopic degrees of freedom is known as the emergent gravity 
paradigm. It has received signihcant amount of support from later investigations, especially 
from the following results: (i) It is possible to express the action functional for gravity as the 
sum of a bulk term and a surface term with a “holographic” relation between them. This result 
holds not only in Einstein gravity but also in all Lanczos-Lovelock theories of gravity [36-38]. 
(ii) More recently, it has been shown [39] that the total Noether charge in a 3-volume TZ related 
to the time evolution vector field can be interpreted as the heat content of the boundary dTZ of 
the volume and the time evolution of the spacetime itself can be described in an elegant manner 
as being driven by the departure from holographic equipartition measured by (Nbuik — Agur)- 
Here, the number of bulk degrees of freedom Abuik is related to the Komar energy density while 
the number of surface degrees of freedom A^sur is related to the geometrical area of the bound¬ 
ary surface. All these results generalize in a non-trivial manner to Lanczos-Lovelock theories of 
gravity [40]. 

These facts suggest that the gravitational field equations on (or near) any null surface in 
any spacetime might have a natural thermodynamic interpretation. One should be able to 
understand all the previous results as different facets of a unified picture and also generalize 
them to an arbitrary null surface. In particular, it should be possible to obtain from Einstein’s 
equations the thermodynamic identity T6\S = 5\E-\-P5\V near any null surface, thereby freeing 
the earlier demonstrations from their restrictive assumptions, like the assumption of spacetime 
being static [15] or having a specific horizon. In this paper, we will provide such a generalization 
of these results for an arbitrary null surface in an arbitrary spacetime, which is neither static 
nor spherically symmetric. 

As we have mentioned earlier, it is possible to attribute thermodynamical entities like tem- 
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perature and entropy to any null surface by introducing local Rindler horizons. Then, with a 
suitably defined pressure P (or more precisely, a work function), we can consider an infinitesimal 
displacement of the horizon along the affine parameter A of the null geodesics off the surface and 
show that Einstein’s equations imply the relation: 

TSxS - PdxV = SxE (1) 

for a suitably defined energy E. Since we already have well-defined, physically motivated, expres¬ 
sions for T, S, P and V, it is now possible to identify the energy associated with the null surface, 
which appears in this thermodynamic identity for an arbitrary null surface under consideration. 
Further, starting from this result and taking suitable limits we can arrive at the previous results 
in the literature as special cases. 

We will obtain the result in Eq. (1) using the component of the Noether current along the 
null geodesics on the surface and the relation between Ricci scalar for the full spacetime with 
the Ricci scalar for the two-surface. This provides a physically well motivated route to study 
thermodynamic structure of the spacetime which has proved to be quiet useful in the past [39,40]. 

The paper is organized as follows: In Section 2, we will introduce the Gaussian null coordi¬ 
nates near an arbitrary null surface which we will use throughout the paper. Then, in Section 3, 
we describe three natural projections of Einstein’s equation which arise in the presence of a 
null surface and their physical interpretation. In Section 4, we introduce the Noether current 
associated with a vector field from a simple geometric identity (which does not require talking 
about diffeomorphism invariance of an action etc.) and identify the vector for which we shall 
compute the Noether current. Next we evaluate the component of the Noether current along the 
null geodesics on the null surface and use it to derive the thermodynamic identity by considering 
a virtual displacement of the null surface along the auxiliary null geodesics. Subsequently, in 
Section 5, we reproduce the previous results available in the literature by specializing to station¬ 
ary and static metrics as well as to spherically symmetric (but not necessarily static) metrics. 
Finally, we conclude with a short discussion of our results. 

Throughout this paper, we use metric signature (—, -f, -f, -f) with the fundamental constants 
G, h and c being set to unity so that Einstein’s equations reduce to Gab = ir^Tab- The Latin 
indices a,b,... run from 0 to 3 and stand for spacetime coordinates, Greek indices ... 

run from 1 to 3 and represent coordinates on the null surface and capitalized Latin indices 
A,B,... stand for coordinates on the 2-surface transverse to the normal to the null surface and 
its auxiliary null vector. 


2 Gaussian Null Coordinates (GNC) 

Since we are interested in the form of the field equations near a null surface, we will begin by 
introducing a coordinate system {u,r,x^) adapted to the null surface. This coordinate system 
will be constructed in close analogy with what we expect in a local Rindler frame and will have 
the following properties: (a) There will be only 6 free functions in the metric thereby freezing 
all redundant gauge degrees of freedom, (b) The null surface we are interested in is chosen to 
be a surface r = 0. Further, r = constant but non-zero surfaces will represent timelike surfaces 
with r —0 leading to the null surface we are interested in as a limit, (c) Observers at rest in 
this spacetime with constant values for {r,x^) will be analogous to local Rindler observers and 
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will perceive r = 0 as their local Rindler horizon. Let us briefly review how such a coordinate 
system can be constructed around any null surface in any spacetime. 

Any arbitrary null surface can be parametrized using Gaussian null coordinates (henceforth 
referred to as GNC), which can be constructed in analogy with standard Gaussian normal coor¬ 
dinates associated with, say, a spacelike surface. In the non-null case, the construction proceeds 
by using geodesics normal to the surface. This construction breaks down in the null case, since 
geodesics with tangent vectors along the surface normal, are actually on the null surface. This 
problem is avoided by introducing an auxiliary null vector satisfying (.ak°‘ = where (.a 
is normal to the surface, and then constructing the coordinates by moving away from the null 
surface along the null geodesics of fc“. The construction of this coordinate system has been 
detailed in [41-43] and we will only recall its essential properties. The line element adapted to 
an arbitrary null surface (identified with r = 0) takes the following form in GNG: 

ds^ = —2radu'^ -\- 2dudr — 2rj3Adudx^ qAsdx^dx^ ( 2 ) 

This line element contains six independent parameters a, Pa and qAB, all dependent on all the 
coordinates (u, r, . The metric on the two-surface (i.e. u = constant and r = constant) is 
represented by qAB ■ The surface r = 0 is the fiducial null surface but surfaces with r = non-zero 
constant are not null. The normal ^a = daX to the r = constant surfaces will be a null vector on 
the r = 0 null surface. Hence, the null normal £“ and the corresponding auxiliary null vector 
have the following components [43] in this coordinate system: 

4 = (0,1,0,0), r = (1,2 to + 7 - 2 / 32 , r/3^) (3a) 

4 = (-1,0,0,0), = (0,-1,0,0) (3b) 

While dealing with an arbitrary null surface with normal there is some freedom in the 
choice of k°'. In our case, we have chosen the auxiliary vector as the vector that was used 
in the construction of the GNG [43] itself, and hence, up to a sign, is a basis vector in the 
GNC. Once the coordinate system adapted to the the null surface is fixed, this particular is 
specified by the conditions = 0 in addition to the usual conditions on an auxiliary vector, 
fc“fea = 0 and = —1- (Here, Ca (with A = 1,2) denote the two coordinate basis vectors 
on the u = constant 2-surfaces on the null surface.) This allows us to work with a physically 
well-defined basis {-^o, 4, e^}- Later on, we will take the projection of certain vectors along 
these basis vectors. The choice made here allows us to take the projections which have direct 
thermodynamic meaning; if we use a linear combination of these vectors or make some other 
choices, then, of course, the projections will get mixed up with each other and one cannot provide 
a simple interpretation to them. 

The non-affinity parameter k for the null normal is defined via the relation = k£°'. It 

turns out that the non-affinity parameter for the null normal we are considering is k = a thereby 
allowing us to interpret the a, which occurs in the metric in Eq. (1), as the surface gravity. The 
vector = —djdr is tangent to the ingoing null geodesic (ingoing since it points in the direction 
of decreasing r), which is affinely parametrized with affine parameter r. We denote Xh to be 
the value of the affine parameter on the null surface. In the remaining discussions, we will work 
with A defined through the following relation: r = X — Xh- It is also useful to introduce the 
vector field: 

C=^ = (I,0,0,0) . (4) 
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which goes to on the null surface. This vector is special since it corresponds to the standard 
time direction in some well-known spacetimes that can be obtained as special cases of GNC (see 
Appendix B) and in the local Rindler frame. We shall describe the dynamics of the null surface 
from the point of view of observers moving along the integral curves of in the region r > 0. 
(We have arbitrarily chosen one side of the null surface. We could have as well chosen the r < 0 
side.) Thus, representing the time direction for our fiducial observers, has to be timelike in 
the r > 0 region. In order to achieve this property for we shall assume that a > 0, at least 
near the null surface in the r > 0 region. (This restriction is also consistent with the idea of 
identifying Q;|r=o to be the surface gravity and associating a temperature aj^'n with the null 
surface.) 


3 Projections of Einstein’s equations 

The vector introduced earlier, when normalized, gives the four-velocity a fundamental 

observer with r,0,<f) = constant in the spacetime described by the metric in Eq. (1). Further, on 
the null surface, —>•£“. Therefore the flux of the matter energy momentum tensor through 
the null surface is determined by the four momentum: 

ga ^ ^ j.a^b ^ ya ^ ^5) 

where the second relation holds in the limit of r —>■ 0. When field equations Gab = SirTab hold, 
we find that 5”“ = = (l/87r)GQb£^ on the null surface. Algebraically, however, it turns out 

to be simpler to concentrate on the Ricci tensor rather than Einstein tensor and define a closely 
related vector: 

= 2Rt£‘’ = 167r2\“£^ (6) 

The structure of Einstein’s equation near a null surface is determined by the components of 
on the null surface. To investigate these components, let us expand in the orthonormal 
basis [43] made of (£“,fc“,e"^) as P°‘ = 4>i£°‘ + (j) 2 k°' + This allows us to construct two 

scalars {(f>i,(j> 2 ) and one transverse vector (pA from the combinations: (pi = —P°‘ka,(p 2 = —P°‘£a 
and the projection P^^q^. (Of these, P°‘ka and P°‘q\ together represent the three components 
of the projection of the flux on the null surface while P°‘(.a brings out the component along 
which is the tangent vector to the ingoing null geodesic.) 

Remarkably enough, these three components {P°'qa, P°'£a, P°‘ka) lead to the three sets of results 
obtained in the literature earlier. The first one {P°‘qa) leads to the Navier-Stokes equation, the 
second {P‘^ia) is related to Raychaudhuri equation and the associated results while the third 
one [P°^ka) leads to the thermodynamic identity in the special cases considered earlier in the 
literature. We will briefly describe these three and then investigate the last one in detail. 

• The contraction P‘^£a' Contraction of the momentum with the null generator of the 
null surface, leads to the standard Raychaudhuri equation, 

Rabti'’ = + a;abCc“'’ - ^, (7) 

involving the combination Rab£°'£^- It is this Raychaudhuri equation which was used in 
the work by Jacobson [32], along with the crucial assumption of vanishing 0 and a for the 
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chosen null surface, to obtain Einstein’s equation from Clausius relation, i.e., 5Q = TdS. In 
other words, the component of responsible for Jacobson’s thermodynamic interpretation 
of Einstein’s equations is obtained by the contraction of P°‘ with the null generator ^a■ But, 
for providing this interpretation, one needs to make additional restrictive assumptions (like 
the vanishing of shear and expansion) which will not hold on an arbitrary null surface. 
Later on, some of these restrictive assumptions in Jacobson’s work were lifted in [44] but 
this demanded the interpretation of the shear and expansion terms as dissipative effects. A 
more detailed discussion of the differences between Jacobson’s approach and our approach, 
as well as some crucial issues in the former are highlighted in [45]. It should also be noted 
that the approach based on the Raychaudhuri equation cannot be generalized in a simple 
manner to Lanczos-Lovelock models, while it turns out that our approach does generalize 
in a straightforward manner to Lanczos-Lovelock theories of gravity [46]. 

We want to work with an arbitrary null surface with non-zero shear and expansion and we 
want to develop an approach which will generalize to Lanczos-Lovelock models in a natural 
fashion. It, therefore, turns out to be more fruitful to study the projection of P“ orthogonal 
to la, especially the contraction P‘^ka- As we shall show, the contraction along has 
a neat thermodynamic interpretation without additional assumptions or introduction of 
dissipative effects. The thermodynamic interpretation of the contraction with P°‘ka also 
generalizes in a straightforward manner to Lanczos-Lovelock theories of gravity [46]. In 
short, the projection on ka leads to richer thermodynamic content. The physical reason for 
this could be the following: Note that, since 1^ = 0 and £°‘ka = —1, the contraction P°‘ka 
actually picks out the components of P“ along £“ which is intrinsic to the null surface. 

• The contraction P°'qa-' Let us start with contraction of P“ with the transverse metric qab 
which is proportional to Rabl°‘q^'^- This expression — when worked out in detail — leads 
to the Navier-Stokes equation on the null surface [33]. More specifically, using vectors and 
derivatives intrinsic to the null surface, the contraction of P“ with qab leads to [33] 

Rmnrq2 = qT£e^m + OTla - Pa + |^ + q^ (8) 

where Da is the covariant derivative defined on the null surface using the projector qab 
and £e denotes the Lie derivative along the null generator 1. We have also separated 
out the trace of Qab and have defined a new object amn = Qmn — It is 

clear from Eq. (8) that it has the form of a Navier-Stokes equation for a fluid with the 
convective derivative replaced by the Lie derivative. This correspondence allows us to give 
the following interpretations to geometric quantities on the null surface: (i) The momentum 
density is given by — fla/87r where Da = nka -f PVjka- In the coordinates adapted to the 
null surface. Da has only transverse components which are given by Da = (3a', this suggests 
interpreting (3a as the transverse fluid velocity. Eurther, we have identified the (ii) pressure 
k/Stt, (hi) shear tensor defined as amn, (iv) shear viscosity coefficient rj = (I/I67r), (v) bulk 
viscosity coefhcient ^ = —I/I67r and finally (vi) an external force Pq = TmaP^- 

• The contraction P°‘ka '■ The contraction with k°', as we have mentioned, has very interesting 
consequences and has not been explored adequately in the literature (except in some special 
cases which we will mention in the sequel). Since this contraction picks out the component 
flowing along the null geodesics on the null surface (i.e the component of P“ along ^“) it 


6 


encodes an intrinsic property of the null surface. It is, therefore, worthwhile to examine 
this in detail for a general case which will be the main thrust of this paper. We will show 
that this leads to the thermodynamic identity we are after. 

The above separation of the components of along provides a clear picture of 

different aspects of gravitational dynamics on a null surface and allows us to identify which of 
the previous results arise from which component of P“. 


4 Thermodynamic identities from gravitational dynamics 
near a null surface 

We are interested in the structure of and its interpretation as a thermodynamic identity. 

To study this, we begin (following [45,47]) by introducing the notion of a transverse metric 
and the work function P. Let Ua be a normalized timelike vector while be another 
normalized but spacelike vector related to our null vectors {la, ka) by Ua = {l/2A)la + Aka and 
ra = {l/2A)la — Aka, where A is an arbitrary function. Then the transverse metric defined 
9ab ~ UaUb — rar^ = lakb + Ibka- The work function of the matter is defined [45,47] as 
P = {l/2)Tabg±‘ = Tabi°‘k^- (In the case of spherically symmetric spacetime, P will be the 
transverse pressure; we will not bother describe the physical meaning of P here since it has 
been done in previous literature.) When Einstein’s equations hold, the work function will be 
proportional to Gabl°‘k^ = {l/2)Gabg°^■ 

We will now study the form of equations which arise when we project the field equations 
along which will lead to a thermodynamic identity. While this can be done directly (see 
Appendix A for such derivation), it is nicer to obtain it from the expression for a Noether current 
which we will now briefly introduce. 

One can associate a natural conserved current with any vector field v°' by 

choosing antisymmetric second rank tensor field corresponding to this vector field as 

167r J“'’ = (9) 

The resulting conserved current J“ is indeed the standard Noether current but this approach 
delinks the Noether current from diffeomorphisin invariance of the action etc. and attributes the 
conservation law to a simple identity in differential geometry. This conserved current J“, for the 
vector field v'^, has the following expression in general relativity [48]: 

I67rJ“(i;) = V6[V“7;'’ - V'’u“] = (10) 

where 

iv“ = -r“ + (1/2) + s^rl) (ii) 

is a linear combination of Christoffel symbols. (Its physical significance is discussed in [48] and 
will not be repeated here). 

For our purpose, we will concentrate on IOtt/co J“(^), which contains the combination Rab^°‘k^ 
(which will become Rabl°‘k^ on the null surface). This is given in our coordinate system by 
— 167rJ“, the component of the Noether current along the null geodesics on the surface. It can 
be worked out in the most general case (presented in Appendix C), but algebraic complexity of 
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the resulting expressions hide the physical interpretation. To bring out the physics involved, we 
will consider a slightly constrained situation in the main text, leaving the discussion of the most 
general case to Appendix C. The simpler case is obtained by setting (a) /3A|r=o = 0 just on the 
null surface but is arbitrary otherwise and (ii) imposing hypersurface orthogonality on the 4- 
velocity constructed out of the vector (see Appendix B). Then, from Eq. (75) in Appendix B, 
we get the result dAQ\r=o = 0 on the null surface. Thus, the two conditions (viz., fiA\r=o = 0 
and hypersurface orthogonality for the 4-velocity constructed out of ^a) lead to the result that a 
is independent of the transverse coordinates on the null surface, which can be thought of as an 
extension of the zeroth law of black hole thermodynamics to a null surface in a time dependent 
situation. 

In this case, the Noether current contracted with ka has the following expression (see Appendix C 
for details): 


IGirkaJ^'i^) = 4:dra + 2adr \n^/q 


( 12 ) 


However, from Eq. (10), we can also rewrite the above contraction of Noether current as follows: 

16^A:,J“(e) = + ka9^^£iKj 

= 2GabCk^ -R + ka9^^£iN- (13) 

We next write the Ricci scalar R in terms of the two-dimensional Ricci scalar R^'^^ for the 
two-surface as 

i (i? - = _2d,^ - 2adr In (v^/3^) - 

+ ^duQABdrq'^^ + du In ^/qdr Iny/q - ‘^dudr^/q ■ (14) 
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Further, the Lie variation term in Eq. (13) has the following expression (see Appendix C): 

ka9^^£iN^j = -2£iK^ - q^^£iNlB = -2dA\n^+ ^^d^qABdrq^^, (15) 

Thus, using Eq. (14) and Eq. (15) in Eq. (13) and using Einstein’s equations along with Eq. (12), 
we obtain the expression for Tab^‘^k^ as 

-TabCk^ = ^( - . (16) 

In the null limit, —Tab^°'k^ = T°'^ka^b = is the work function defined previously [47], which 
enables us to write the above equation on being multiplied by ^/q as 



Stt 


dxduy/q 


(17) 


Multiplying the above equation by 6X = 6r and then integrating over a, u = constant slice of the 
null surface with area element d^Xy/q, we arrive at 


J d^Xy/^SXT^^ = 





-Sx 


',2 



(18) 
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where we have made use of dAOt\r=o to take a outside the integral and used the notation 5\ = 
{8X)d\. {5\ can be thought of as the change due to a virtual displacement along the vector 
= —k‘^SX. We shall explain the meaning of virtual displacement shortly.) 

We take (a/27r) as the temperature of the null surface related to the surface gravity a and 
dS = d^Xy/q/A as the entropy associated with a proper transverse area element d^Xyfq. Further, 
in the relation 


J Py/qd^xSX 


j d^Xy/^SX = J d^Xy/^SX T'" 


(19) 


where P is the work function, the integral on the left hand side represents the amount of work 
done in a virtual displacement of the null surface by an amount 6X along ka- Then Eq. (18) can 
be recast in the following form; 


FSX = TSxS - 6xE , (20) 

where we have defined the energy swept out by the null surface, S\E, as 

S,E = M (I - (21) 

with X representing the Euler characteristic of the 2-dimensional slice of the null surface trans¬ 
verse to £°' and 


X=^l ( 22 ) 

and A± representing the area of this slice. (If the 2-dimensional surface is not compact, we cannot 
introduce the Euler number but x is still defined by this integral in our result.) Performing an 
indefinite integral along A, this can also be written as: 

As an aside, we note that when /3^ ^ 0 we pick up a ‘kinetic energy’ term {1/2)(3aP^ in the 
expression for the energy and the result is given by (see Eq. (54) in Appendix A) 

The notion of the virtual displacement introduced here is a straightforward generalization of 
the idea discussed in previous works (see e.g. [7,15]). (In these earlier works, one considered 
spherically symmetric and static spacetimes and concentrated on the horizon as the null surface; 
here we have made no restrictive assumptions and deal with an arbitrary null surface.) To see the 
correspondence explicitly, consider the simpler situation of a static and spherically symmetric 
spacetime, such that —gtt = = /(*")j with a (non-extremal) horizon at r = a. Then /(a) = 0 

with f'{a) related to the horizon temperature by T = /'(a)/47r with ‘prime’ denoting derivative 
with respect to r. Repeating our exercise, treating the horizon as our chosen null surface, will lead 
to the relation: f'{a)a — 1 = SttPo^, where P = T^T is the radial pressure. (The analogue of this 
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relation in the general case of an arbitrary null surface is given by Eq. (17)). If we multiply this 
relation by 5a, we can rewrite the equation in the form T5S — SE = P6V, purely algebraically. 
We can interpret this relation in this case by considering two solutions of the field equations 
differing infinitesimally such that horizons are located at a and a + Sa with all other infinitesimal 
differences treated as the differences between these two solutions. (This is analogous to the 
relation for general null surface given in Eq. (20)). Hence, the virtual displacement is essentially 
a shift between the location of the fiducial null surface in two solutions of field equations. The 
shift moves the null surface by an amount (5A, where A is the affine parameter along the null 
geodesics of /c“. More detailed discussion of this idea can be found in [9,12,15,26]. 

Coming back to the case with Pa = 0, we note that F represents the integral of the work 
function over the null surface i.e. F = f P^dPx. Then Eq. (20) is better interpreted when 
rewritten in the form 


5xE = TSxS - F5X . (25) 

This expression is quite suggestive. The virtual displacement can be interpreted as a physical 
process that displaced the null surface from r = 0tor = r + 6X. Then, the energy engulfed 
by the null surface in this displacement is the sum of a heat energy (viz. the temperature 
multiplied by the change in entropy) and the work done during this virtual displacement of the 
null surface. The above equation can also be interpreted as the total energy in the region being 
a sum of an energy corresponding to matter, represented by the work done term, and energy 
corresponding to pure gravity, represented by the heat term. This is the most general form of 
the thermodynamic identity which arises from the projection of along k°'. We shall now 

discuss the applications of this result to special cases. 


5 Special Cases 

In the previous section, we have shown the equivalence of gravitational Held equations with a 
thermodynamic identity for an arbitrary null surface. The result in Section 4 has been obtained 
in the earlier literature for some special cases. In this section, we will connect up with the 
earlier work by specializing this result to (i) stationary spacetimes without any other symmetry 
(ii) spherically symmetric spacetime which is not necessarily static and (iii) static spherically 
symmetric spacetime. 


5.1 Stationary spacetime 


Since we have identified as our time flow vector, stationarity involves setting partial derivatives 
of metric components with respect to u to zero (see Appendix B). In this case the thermodynamic 
identity in Eq. (20) retains its form with a simpler expression for the energy term. The expression 
for energy in Eq. (21) becomes 


5xE = = ^5r 

4 4 

This immediately implies 

dE _x 
dr 4 


(26) 


(27) 
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This matches with the result in [15]. Notice that, even in the more general case (when E is 
given by Eq. (24)), we can obtain the same result if we assume (a)/3^ = 0 on the null surface 
and (b) the stationarity condition, viz. the metric is independent of u. So the above result does 
not require the spacetime to be static (which involves the additional condition of hypersurface 
orthogonality) but only requires stationarity. 

The additional restrictions required for achieving staticity are the conditions of hypersurface- 
orthogonality given in Eq. (73) and Eq. (74) in Appendix B, reproduced below: 

r (/329./3 i - PidrP 2 ) - [dih - 92/3i) = 0 , (28) 

radr^A — dACt — = 0 . (29) 

As we have noted in Appendix B, these conditions would imply dAOi\r=o = 0. There is no 
modification to our result above since we had already assumed hypersurface-orthogonality. 


5.2 Spherically symmetric spacetime 

To restrict the GNC metric in Eq. (2) to a spherically symmetric form, the most convenient 
way would be to enforce the geometry of 2-spheres on the u =constant, r=constant 2-surfaces. 
However, the u =constant surfaces should not be considered as constant-time surfaces as these 
surfaces are actually null. Thus, identifying with the angular coordinates, we should demand 
dACt = 0, = 0 and qab = f{u,r)dfl‘^ = f{u,r) {dO"^ + s\t? Odcf)^) to arrive at a spherically 

symmetric form of GNC. The form of the line element now becomes 

ds^ = —2ra{r,u)dv? + 2dudr + f(u^r)dVt^, (30) 


which is of spherically symmetric form [49] . The 2-surface we are mainly interested is u =constant 
cross-section of the null surface at r = 0. Defining the “radial coordinate” [50], R{u,r) = 
\Jf{u, r), and expanding in Taylor series in r around r = 0, we obtain i?(u, r) = RH{u)+rg{r, u), 
where the last term is not just the linear order term in the Taylor expansion but represents all 
the higher order terms taken together. The null surface at a constant u has a “radius” Rh{u), 
with the M-dependence allowing for the area of the 2-surface to be changing with u. (So we have 
assumed spherical symmetry but have allowed for time dependence). 

Again, our result holds with a simpler expression for the energy. To see this, let us look at 
Eq. (21). For 2-spheres, the Euler characteristic equals 2. Substituting in Eq. (21), we obtain 


S\E = SX 


Stt J ’ 


(31) 


where A± is the area of the compact 2-surface. Interpreting E as the energy associated with a 
u =constant, r=constant 2-surface, we have E = E{u, r) for the spherically symmetric spacetime. 
The above equation can then be written as a partial differential equation as follows: 

dE _l dxduAx 

aA “ 2 Stt ’ ^ ^ 


where d/d\ is the same as djdr, since r = X — Xh, and is taken keeping u constant. The solution 
is 


E{X, u) 


X 

2 


duA± 

Stt 


E{u), 


(33) 
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where F is an arbitrary function. In order to fix this function, let us consider the ingoing null 
geodesics, —djdr on the u =constant surface. These are 9 =constant, 4> =constant lines and 
are hence radial. Let us assume that moving along this geodesic will lead us to intersect ever 
smaller two-spheres, i.e we will be moving towards the inner part of the 2-spheres. (If this is not 
the case, we can move in the other direction.) Due to spherical symmetry, all these geodesics 
will intersect at the common center of the 2-spheres (assuming that the geometry is such that 
this center exists). Since all the geodesics are given the affine parameter value Xh at the r = 0 
2-sphere, they will all have the same value of affine parameter A (or, equivalently, r), at the 
center. But this affine parameter value may be different on a different u-constant surface. Let 
us label this value Xo{u). At the center, Eq. (33) becomes 


E{Xo{u),u) 


Xo{u) 

2 


duA± 

Stt 


X=Xo{u) 


F{u), 


(34) 


Now, A± = Fnj = AttR^ so that = SnRduR, which is zero at A = Ao(m) since R = 0. 

Thus, we obtain 

E{Xo{u),u) = + F{u). (35) 

Since it seems natural to associate zero energy to a 2-sphere that is essentially a point, we shall 
choose F{u) = — Ao(w)/2. Thus, our definition for energy becomes 

T?f\ \ X-Xo{u) duA± , . 

E{X,u) = - (36) 

For the r = 0 surface, we have the energy as 

,, , Xh-Xo{u) duA±\^^^ 

EiXn.u)- -^-- 

= - R„a,R„ . ( 37 ) 


A special case of interest is the case in which the 2-sphere line element is (r -|- Rh)^(EX'^, for a 
constant Rh, which puts i? = r -|- Rh- Taking i? as a coordinate instead of r, we obtain the 
metric element 

= -2{R - RH)adu^ -f 2dudR -f R'^dfl^, (38) 

which is of the form of the metric element used in [51]. (The metric that the authors of [51] start 
with is different in form, but they reduce it to this form through a coordinate transformation.) 
In this case, duRn = 0, Xh — Xq{u) = th — rQ{u) = Rh — Rq = Rh, where we have used 
r = A — Xh, denoted r at our fiducial null surface as rn, even though it is zero in our framework, 
and the r and R at the center of the 2-spheres as tq and i?o respectively. Substituting, we get 
the energy of a u =constant, r = 0 surface as 

E{Xh,u) = ^ (39) 


matching with the results of [51]. However, we should stress that they have derived the first law 
for an apparent horizon which is not a null surface while in this work we have derived a similar 
result for a null surface. 
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With hindsight, it may have been better to set the origin of the afhne parameter at the center 
of the 2-spheres, i.e Ao(u), to obtain a form for the energy as 

E{\h, u) = - RhOuRh , (40) 

where A//, the value of the affine parameter at r = 0, is no longer a constant but depends on u. 


5.2.1 Stationary spherically symmetric case 

Since = 0 and 9^0 = 0, enforcing the stationary conditions d^a = 0 and d^f = 0 leads 
to staticity (see Eq. (28) and Eq. (29)). This is not surprising as every stationary spherically 
symmetric spacetime is automatically static [50]. Thus, the previous results for static spacetime 
hold with the spherically symmetric transverse metric. In this case Eq. (37) leads to 


while using Eq. (40) leads to 


E = 


Ag — Ao(u) 
2 


E{Xh,u) 


Xh 

~Y' 


(41) 


(42) 


which matches with previous results [15]. Note that the result in [15] was obtained by setting 
an arbitrary integration constant to zero to ensure E —^ 0 when Xh = 0. In the spherically 
symmetric case, this has a very physical interpretation as the radius of the 2-sphere shrinks to 
zero. 


6 Discussion 

We started this work trying to address the question of whether the gravitational field equations 
near any null surface in an arbitrary spacetime reduces to a thermodynamic identity, generalizing 
results previously available in the literature for special cases. We have shown in Section 4 that 
this is indeed possible, by introducing (a) the temperature through surface gravity, (b) entropy 
density from the area and (c) the work function from the transverse metric We then 

obtain, by projecting the Einstein’s equations along the direction, a relation of the form 
T5\S = 5xE + P5\V where the variations represent virtual displacements of the null surface 
along null geodesics off the surface. 

Given an arbitrary null surface with associated normal (£“), co null vector (fc“) and the 
transverse metric (g^), one can study the projections of the vector P“ oc R^t ^ along each of 
these. We pointed out that the projection of along g[[ leads to the Navier-Stokes equation 
on the null surface while the projection along is related to the Raychaudhuri equation. This 
clearly shows that all the information contained in the field equations posses thermodynamic 
interpretation. 

As an aside, note that our result, arising from projection P°'ka along fc“, is distinct from 
any result (like e.g the connection with Clausius relation [32]) obtained from projection P°‘ia 
along (and the resulting the Raychaudhuri equation) and these two class of results should 
not be confused with each other. As noted earlier, the results based on projection along 
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(viz. the Raychaudhuri equation) could not be generalized in a simple manner to Lanczos- 
Lovelock models and they need additional restrictive assumptions (like the vanishing of shear 
and expansion) even for thermodynamic interpretation. In fact, as the earlier work in [45] clearly 
points out, the thermodynamic structure of the curvature tensor is not properly captured in the 
components which occur in the projection along and the Raychaudhuri equation. We believe 
the other two projections (on qab and on ka) leads to richer thermodynamic content. This is 
because they pick out the components of along (due to = 0,£‘^ka = —1!) and along 
both of which are intrinsic to the null surface. They are also most likely to remain valid even 
in a more general class of theories. (We already know that if the spacetime is static, then the 
resulting thermodynamic identity holds even for Lanczos-Lovelock models.). This clarification 
of the different projections of and their thermodynamic relevance is an important offshoot 
of our work. 

We derived our result starting from the Noether current, which shows again the intimate 
connection between Noether charge and thermodynamics seen in earlier works. Through this 
exercise, we have introduced a definition of energy which reduces to energy definitions introduced 
previously in the static case [15]. In the most general context, this involves time derivatives of 
the area of the null surface and additional terms involving off-diagonal metric elements /3 a- 
If we assume /3a = 0 on the null surface, enforcing hypersurface-orthogonality on our chosen 
time-flow vector naturally leads to dAT\r=o = 0, which is an extension of the zeroth law of 
thermodynamics to the case of an arbitrary null surface. In this situation, the energy consists 
of two parts: the standard two dimensional curvature scalar, related to the Euler characteristic 
of the null surface, and a term involving time rate of change of the null surface area. Since the 
two-metric is independent of time in the static case, energy becomes solely dependent on the 
Euler characteristic of the null surface. We then discuss the case of stationary, static, spherically 
symmetric and stationary spherically symmetric spacetimes and make connection with results 
previously available in the literature. 

To summarize, we have shown that for any arbitrary spacetime, without assuming any sym¬ 
metry, gravitational field equations in general relativity near an arbitrary null surface reduces to 
a thermodynamic identity. Also, for a restricted class of spacetimes with hypersurface orthogo¬ 
nality enforced, zeroth law holds (even in time dependent cases). It is interesting to ask whether 
identical results hold for the most general class of gravitational Lagrangians with second order 
equations of motion, i.e. in Lanczos-Lovelock gravity, as well since the previous result for static 
spacetime was indeed applicable to Lanczos-Lovelock models. This work is under progress and 
the results will be presented elsewhere. 
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Appendices 


A Derivation from gravitational field equations 

In a static and spherically symmetric spacetimes, we have the relation G\ = G^. between the 
Einstein tensor components in the near horizon limit. This relation continues to hold for arbitrary 
static spacetimes as well [52]. Taking a cue from this, let us evaluate the corresponding Einstein 
tensor components for the GNC coordinates. All expressions are evaluated on the null surface, 
i.e at r = 0. We have 

u — y ^au — ^UT'i — y ^ ar — ^ur 

Gl = Gl= Gur = Rur -lR = Rur-l {2Rur + Rab) = RaB 

= ^aq'^^drQAB - + ^q^^drduqAB + dr\n^duln^+ ]^duqABdrq^^ 

+ ^ [13^Pa + P^q^^dAqcD - 2P^q^^dcqAD + 2q^^dAPB) 

= adr In ^ + ^q^^drduqAB + dr In y/qdu In + ^duqABdrq'^^ 

+ i {P^Pa + P^q^^dAqCD - 2p^q^°dcqAD + 2q^^dApB) (43) 

Also, we have the following identity: 

-\q^^RAB = {Rll + RTd5E + Rll) 

= —SgR^'^ — -ScdRab^ (44) 

where = SqS^ — S^Sq and we have used the relation 

SsRTr = SiRuAr = q^"^ RbuAt 

= q^^RBuAr = q^^RArBu = S^R^ (45) 


The terms in Eq. (43) involving qAB and its derivative can be manipulated leading to 


-q drduqAB + dr In y/qdu In + -duqABdrq 


AB 


= (q^^^duqAB) - \9rq^^ 
_ (du\/q\ I dryfqduy/q 


duqAB + dr In y^du In + -duqABdrq^^ 
- "^drdu-^^ 


(46) 
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The terms with (3a in Eq. (43) can be simplified leading to 

i {j3^I3a + (S^q^^dAqcD - 2(3^q^^dcqAD + 2q^^dAl3B^ = In V? 

- \dcl3^ + ^qAodc + ^q^^dAiSs 

Then, substituting all the expressions in Eq. (43), we arrive at the following result: 

/^r _ sA jduB ^ xAB jjCD _ ^ ^AB p 

Er- — —Ob^uA ~ '^^CD^AB ~ ~ 2^ ^AB 

= adr In + -^drdu^/q + (48) 

For displacement of the null surface by an amount 5\ along the ingoing null geodesic, we multiply 
the above equation by SX^/q and use the gravitational field equation G(. = SttT^T, leading to 


8TTT((6X^/q = aS\^/q 


-R^^^SXy/q + SXdxduy/q - 


+ 2 ^ 3 . (AiP^) 


SXy/q 


(49) 


where we have used the relation, S\f = {df/dX)6X, for any scalar function /. Then, dividing 
the above equation by Stt and integrating over a two dimensional surface with d?x we arrive at 


j d^x^dxr; = I d^x (£) (5a - 5X{ I d^x^^^ I d^xdxdu^ 


d^ 




-/3^ + — 

^ 2yg 


dAiVdP^) } 


(50) 


The null rays tangent to the null surface have the non-affinity coefficient a, which suggests 
defining (Q:/27r) as the temperature of the null surface. Along with this, we can interpret T(( as 
the normal pressure P± on the null surface. This identification allows us to interpret the object 


F = 


/ 


d?Xy/qPx 


(51) 


as the average normal force over the null surface. Then, FdX can be interpreted as the virtual 
work done in displacing the null surface by (5A along ingoing null geodesics. Eq. (50) can now be 
written as 

F6X = J d^xTdxs-6xE, (52) 

where s is the entropy density of the null surface with the following expression: s = (y/q/4), 
which equals the Bekenstein-Hawking entropy density. We have also identified the energy E 
associated with the null surface as 

E=^ldxl I |%A I + 

(53) 
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When the two-dimensional surface is compact, this reduces to a simpler form, given by 

where x represents the Euler characteristic of a two-dimensional compact manifold M 2 without 
boundary and is given by the following expression: 

x{M2) = ^[ (55) 

M Jm2 

Note that in this most general situation the first law has to be interpreted as follows: under 
infinitesimal shift of the null surface along ingoing null geodesics, change in energy and work 
done due to pressure adds up and yield T5xs integrated over the null surface. Thus, in this 
general case TSxs has to be interpreted locally as being due to displacement of a small element 
on the null surface. This difficulty arises since the temperature a/27r in the Tds term is dependent 
on the transverse coordinates and cannot be taken outside the transverse integral. The above 
discussion outlines the derivation of first law from the field equation perspective and matches 
exactly with the one obtained from Noether current formalism in Section 2. 


B GNC metric in static form 

The GNC line element is 

= —2radu^ + 2drdu — 2rfiAdudx^ + qAsdx^dx^ . (56) 

We shall attempt to reduce Eq. (56) to the form of the static metric in [52], 

ds^ = —N'^di^ -\- dri^ + aAndy^dy^ , (57) 

using appropriate restrictions and coordinate transformations. We shall place the first restriction 
on a, demanding it to be positive in the region r > 0. The utility of this restriction will be clear 
in due course. 

Eq. (56) represents the line element near an arbitrary null surface in an arbitrary spacetime. 

To get to Eq. (57), we need to enforce staticity. A static spacetime should satisfy the following 

two requirements [50]: 

i) There must exist a timelike vector that satisfies the Killing condition, i.e. 

Va6 + Vh^a = 0 . (58) 

ii) must be hypersurface-orthogonal. By Frobenius theorem, this is equivalent to demanding 

^[aVa6] = 0. (59) 

If only the first condition holds, then the spacetime is called stationary. For stationarity, it is not 
necessary that the vector be timelike everywhere in the spacetime. If we impose that be 
timelike everywhere, then even the Schwarzschild spacetime with the Killing vector = d/dt 
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will not be stationary. Thus, we will only demand that be timelike in possibly only in part of 
the spacetime (see [53]). 

Since the line element Eq. (56) has been constructed in a region near the null surface and 
no claim has been made about its validity for the entire spacetime, we shall adapt the above 
criteria to our situation by calling a GNC metric as static if we can find a timelike vector that 
satisfies Eq. (58) and Eq. (59) in the region of validity of Eq. (56). We shall further restrict the 
domain of validity to the r > 0 region, where guu < 0 for a > 0, since even in a Schwarzschild 
spacetime the timelike Killing vector is timelike only outside the horizon. 

The next logical step would be to choose a timelike vector in the chosen domain and 
demand that it satisfies Eq. (58) and Eq. (59). While these two conditions are enough to render 
the spacetime static, the static line element in Eq. (57) also has a Killing horizon at n = 0. In 
other words, the norm of the Killing vector vanishes at n = 0. We would like our null surface at 
r = 0 to go to the Killing horizon in the static limit. The Killing vector for Eq. (57) lies on the 
Killing horizon. Thus, we are looking for a vector that is timelike in the region r > 0, is null 
at r = 0 and lies on the null surface r = 0. An obvious choice is the vector ^ = d/du. 

To strengthen the motivation for this choice, we shall now demonstrate that it corresponds 
to the timelike Killing vector in Schwarzschild and Rindler metrics. Both Schwarzschild and 
Rindler metrics have the form of the /(r)-metric: 

ds^ = -f(r)df + + QAsdx^dx^ , (60) 

/(O 

with /(r) = 1 — 2M/r giving Schwarzschild and /(r) = —2Kr giving Rindler. The timelike 
Killing vector in the coordinate order {t,r,x^,x'^) is = (1,0,0, 0). Defining a new coordinate 
u by the relation 

f dr 

we have 

(IT (IT 

du = dt+ dt = du- —— . (62) 

f{r) f{r) 

Substituting in Eq. (60), we obtain 

ds^ = —f{r)du^ + 2dudr + qAsdx^dx^ , (63) 

which is the GNC line element Eq. (56) with /3^ = 0. In coordinates (m, r,x^), = 

(1, 0, 0,0) = d/du^ our chosen timelike Killing vector for GNC. 

Having chosen a we shall now apply the conditions Eq. (58) and Eq. (59). The Killing 
condition gives 

£^9ab = 0 => CdcQab + gcbda^'^ + gacdbC = 0 ^‘^dcgab = = 0 . (64) 

Thus, the Killing condition demands that all the metric components be independent of u. This 
means 

dua = 0; OuPa = 0; duqAB = 0 . (65) 

Next, let us look at Eq. (59). The equation ^[a^a^b] = 0 gives four equations corresponding 
to {a,b,c) = {u,x^,x‘^), {a,b,c) = {r,x^,x^) and two equations with {a,b,c) = {u,r,x^) for 
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A = 1,2. These correspond, respectively, to 


2r^a (di^2 - d2^i) + (/329„,Si - PiduP2) + 2r^ WAa - P2dia) = 0 , (66) 

{hdrPi - jiidrh) - r {di(32 - 92/?i) = 0 , (67) 

2r'^adrPA — 2rdAa + tBuPa — 2r'^PAdrO' = 0 . (68) 

Note that these equations are not all independent. For example, enforcing the last two equations, 
the first one can be seen to be satisfied identically. Thus, we just need to demand 

{hdrPi - Pidrh) - r {dih - 92/3i) = 0 , (69) 

2r^adrPA — 2rdAQ. + rdul^A — 2r^/3^9r.a = 0 . (70) 


These equations are automatically satisfied at r = 0. Elsewhere in the spacetime region under 
consideration, we can cancel a factor of a power of r to get conditions on the metric components. 
Since we are considering smooth functions, we should expect these conditions to hold even at 
r = 0. Thus, canceling off overall factors of constants and powers of r, we obtain the following 
conditions for hypersurface-orthogonality: 


r (/329r/3i - /3i9r/32) - (9i/32 - 92/?i) = 0 , 

2radrl3A — 2dAa + OuPa — 2r/3^9rQ: = 0 . 

(71) 

(72) 

Specializing to the stationary case, we enforce Eq. (65) and obtain 


r {/32drl3i - I3idrl32) - (9i/32 - 02 ^ 1 ) = 0 , 
radr^A — dAOi — rf^AdrCt = 0 . 

(73) 

(74) 

In particular, Eq. (74) implies 


II 

0 

II 

0 

(75) 

Thus, imposing staticity on the GNC metric with ^ = d/du, we get a generalization of the zeroth 
law of black hole thermodynamics. 

Once staticity is imposed, it is advantageous to transform to a coordinate where it is manifest. 
Let us take the hypersurfaces to which is orthogonal to be level surfaces of a function t, i.e 

we shall take to be orthogonal to t = constant surfaces. Then, we should have 

= F(u,r,x^)Vat ■ 

(76) 

We shall show that there exists a t which satisfies this equation if we take F{u,r,x^) = —2ra. 
With this choice of F, Eq. (76) becomes gau = —2raS/at, where Vat = —{gau/‘2ra). Hence the 
components of the vector Vat in GNC coordinate reads. 

{dut,drt,dAt) = {-2ra, 1, -t^a) 

2ra 

(77) 

which immediately leads to an expression for dt as: 


dr PAdx^ 

dt = du -1-. 

2ra 2a 

(78) 
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For dt to be a perfect differential, the following integrability conditions need to be satisfied: 


dA 




Ba 




( 79 ) 


It can be verified that these integrability conditions are satished courtesy the hypersurface- 
orthogonality conditions, Eq. (73) and Eq. (74), that we have imposed. If we transform to 
coordinates we will have 

= duf . (80) 

Hence, the stationarity condition Eq. (65) in coordinates (t, r, becomes 

dta = 0; dt^A = 0; dtqAB = 0 . (81) 


dtf\r,xA = duf 


du 


We are now ready to write down the line element in the static coordinate system. From Eq. (78), 
we have 


du = dt + 


dr PAdx^ 


2ra 2a 

which when substituted in Eq. (56) gives the line element 


(82) 


= —2radt^ + — - ^^drdx^ -I- 

2ra a 


QAB 


dx^dx^, 


2a J 


(83) 


with dta = 0, dtfiA = 0 and dtQAB = 0. This is the GNC line element written in an explicitly 
static coordinate system. 

To transform to the static coordinate system of [52], we can first identify = 2rQ;. The 
second step would be to install a Gaussian normal coordinate system in the spatial slice by 
sending out normal geodesics from the r = 0 surface. (The explicit coordinate transformations 
to reach this coordinate system, however, is difficult to obtain in closed form.) 


C Derivations of expressions used in text 

Let us evaluate a couple of expressions we require in Section 4. Noether current for has the 
following expression J“(^) = Vf, [J“^ {^)] with (^) = We shall make use of the 

following expression from [48] : 

= -v.v,r + ^ • (84) 

The object 16nkaJ°'{£,) can be evaluated most easily by using the following identity for any two 
vector fields and v°‘: 

167™oJ“(z;) - Uag'’^ £vN°-j = 2RabU°-V^ = 16TTVaJ°‘{u) - Vag"^£uNtj (85) 
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Applying the above result for the vectors and we obtain 


= IBTTfc, J“(e) - kag^^£^N- (86) 

= l&i:^ar{k)-ia9^^£kNtj (87) 

We have chosen the auxiliary vector fc“ such that = — VqU. Thus, the Noether current for fc“ 
vanishes. Hence, the above equation can be written to yield the value for 167rfca J“(^) as 

J“(0 = kag^^£iN- - iag^^£kN^^ (88) 

Using Eq. (84), the two Lie variation terms can be calculated in a straightforward manner leading 
to 

k,g^^£^N- = -2£^N:^ - g^^£^NlB = -2dA In v^+ \duqABdrq^^, (89) 

^ag'^£kN^j = g''^ £kNlj = - 2adr In ^ - 2dudr In + ^duqABdrq'^^ 

-/3" - . (90) 


Using these expressions, we arrive at 


167rfco J“(C) = 4:dra + 2adr In Vg + H- —Ba iJqP^) ■ 


(91) 


Now, the Noether current expression can be written as 

2RabCk^ = 167rA:, J“(e) - kag^^£^N- (92) 

The above relation can also be verified by calculating directly in GNC coordinates: 

_p _ r? _ f) P® _ f) P*^ -L P^ P^ _P“ P^ 

^ -£^ur ur ar ur^ ab ^ ub^ ra 

= drK, + dAvt - dudr In \n ^ - T^.T^, - rc.r:;, - 

= -2dra - i/3^ - ^9aI 3^ - dudr In ^ - adr In ^ - ^^P^dA In y/q 

_ -plA _ P-^ P^ _ P"^ P-® 

^uA^ru ^ ur^rA ^uB^rA 

= -2dra - \^dA^^ - dudr In ^ - adr In - ^Id^dA In + i^^gAs^r-g^^ - i/3^ 

(93) 

Expression for Rabi°‘k^ can be rewritten by using the result Rab^°‘k£ = Gab^°‘k£ — {1/2)R and 
the field equation Gab = ^irTab as 

-TabCk^ = -jL [leirkariO - ka9^^£iN-] - A (94) 
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Since energy was defined in terms of Ricci scalar on the 2-surface for arbitrary static spacetimes 
[15], we need a relation between R and Using which we arrive at 



^ - 2dra - - dudr \n^- adr In In -h ^dugABOrq"^^ 

OTT \ Z z 4 

+ 2dra + 2adr \n^+ ^5^ (^/3^) + 





(95) 


Then, integrating over the u = constant, r = constant two dimensional surface with integration 
measure (Pxy/q, we arrive at 



(96) 


where the energy is defined as 


I d^xdr^(^ - 



^^dA (yg/3^) - (97) 


which exactly coincides with the result derived in Appendix A. 
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